The last stand before MAP: cosmological parameters from lensing, CMB and galaxy 

clustering 
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Cosmic shear measurements have now improved to the point where they deserve to be treated on par 
with CMB and galaxy clustering data for cosmological parameter analysis, using the full measured 
aperture mass variance curve rather than a mere phenomenological parametrization thereof. We 
perform a detailed 9-parameter analysis of recent lensing (RCS), CMB (up to Archeops) and galaxy 
clustering (2dF) data, both separately and jointly. CMB and 2dF data are consistent with a simple 
flat adiabatic scale-invariant model with Qa = 0.72 ± 0.09, /i 2 f2 c dm = 0.115 ± 0.013, h 2 Q.b — 
0.024±0.003, and a hint of reionization around z ~ 8. Lensing helps further tighten these constraints, 
but reveals tension regarding the power spectrum normalization: including the RCS survey results 
raises as significantly and forces other parameters to uncomfortable values. Indeed, as is emerging 
as the currently most controversial cosmological parameter, and we discuss possible resolutions of 
this as problem. We also comment on the disturbing fact that many recent analyses (including 
this one) obtain error bars smaller than the Fisher matrix bound. We produce a CMB power 
spectrum combining all existing experiments, and using it for a "MAP versus world" comparison next 
month will provide a powerful test of how realistic the error estimates have been in the cosmology 
community. 



I. INTRODUCTION 

An avalanche of precision data has revolutionized our 
ability to constrain cosmological models and their free 
parameters in recent years [1-36]. Around 1998, cosmic 
microwave background (CMB) data had become good 
enough to allow a realistic high-dimensional parameter 
space to be explored by fitting theoretical models di- 
rectly to the measured CMB power spectrum {e.g., [8,7]). 
In contrast, other cosmological constraints were included 
merely as priors on individual parameters, for instance 
the baryon density from Big Bang nucleosynthesis, the 
matter density from galaxy cluster abundance and the 
Hubble parameter from direct observations. 

With the advent of improved galaxy redshift surveys 
like PSCz [37], 2dF [38] and SDSS [39], the galaxy power 
spectrum P(k) was upgraded and admitted to the "CMB 
club" : it was so accurately measured that people started 
fitting models directly to the measured P(k) curve rather 
than merely to some phenomenological parametrization 
thereof (say in terms of an amplitude and a "shape pa- 
rameter") [8,9,23]. The Lyman a forest (LyaF) power 
spectrum has now undergone the same upgrade [11,30] 

One promising cosmological probe is still conspicuously 
absent from this club and has so far been left out in the 
cold: the cosmic shear power spectrum measured with 
weak gravitational lensing. Like the CMB, it has the ad- 
vantage of probing the dark matter distribution directly 
[40], without murky bias issues. Despite the spectacu- 
lar progress since this elusive signal was first detected 
in 2000 [41-44], increasing the sky area covered by al- 
most two orders of magnitude [40,45-52], it has so far 
only been included in cosmological parameter studies as a 
prior on two parameters: the power spectrum normaliza- 
tion and the matter density. Yet it has been shown that 



the detailed shape of the weak lensing power spectrum 
P K (£) potentially contains almost as much cosmological 
information as the CMB [53-55] , with great prospects for 
degeneracy breaking and cross-checking. 

The goal of this paper is to treat weak lensing on an 
equal footing with CMB and galaxy clustering, fitting 
theoretical models directly to the CMB, LSS and lensing 
power spectra. Our results can also be viewed as the last 
stand before MAP: using our constraints combining the 
state-of-the-art available data for a "MAP versus world" 
comparison next month will provide a powerful test of 
how realistic the error estimates and underlying assump- 
tions have been in the cosmology community. 

The rest of this paper is organized as follows. We 
present our analysis method in Section II, describe the 
CMB, galaxy and lensing data used in Section III, com- 
pute quantitative model constraints in Section IV and 
discuss our conclusions in Section V. 



II. METHOD 

In this section, we briefly describe our method for con- 
straining cosmological models. 



A. General method 

We use the grid-based multiparameter analysis method 
described in detail in [8,9,30], with modifications as de- 
scribed below. It involves the following steps: 

1. Compute CMB and galaxy power spectra and the 
weak lensing aperture mass variance for a grid of 
models in our 9-dimensional parameter space. 
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2. Compute a likelihood for each model that quantifies 
how well it fits the data. 

3. Perform 9-dimensional interpolation and marginal- 
ization to obtain constraints on individual param- 
eters and parameter pairs. 

As in [30], the we parametrize cosmology with the 11 
parameters 



P = (t, fik,fiA,Wdm,Wb,/i/,n s ,nt, A s ,r, b). 



(1) 



They are the reionization optical depth r, the primordial 
amplitudes A s , rA s and tilts n s , n t of scalar and tensor 
fluctuations, the bias parameter b defined as the ratio 
between rms galaxy fluctuations and rms matter fluctu- 
ations on large scales, and five parameters specifying the 
cosmic matter budget, curvature f^k, vacuum energy Ct\, 
cold dark matter f2 c d m , hot dark matter (neutrinos) Q, v 
and baryons Ob- The quantities uJb = h 2 ilb and uj^m = 
/i 2 Odm correspond to the physical densities of baryons 
and total (cold + hot) dark matter (fidm = ^cdm + &v), 
and /„ = fi^/fidm is the fraction of the dark matter that 
is hot. We assume that the galaxy bias b is constant on 
large scales [56] . We make A s as a discretized parameter 
to marginalize and find a best fit of as later in compu- 
tation, as comparing to an undiscretized A s in previous 
papers [30]. Throughout this paper, we assume a neg- 
ligible contribution from spatial curvature and massive 
neutrinos, setting £lk = f v = 0. Our final parameter 
grids are as follows: 

• t = 0,0.05,0.1,0.2,0.3,0.5,0.8 

• O a = 0,0.1, ....,1.0 

• uj dm = .02, .05, .08, .11, .13, .16, .20, .50 

• u h = .003, .015, .018, .020, .022, .025, .03, .04, .07 

• n 3 = 0.5,0.7,0.8,0.9,1.0,1.1,1.2,1.4,1.7 

• th = -1.0, -0.7, -0.4, -0.2, -0.1, 

• r = 0, 0.1, 0.2, 0.3, 0.4, 0.6, 0.8, 1.0, 1.4, 1.8 

• A s = 0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9 

• & = 0.5,0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4 

For each grid point in parameter space, we compute 
CMB and matter power spectra using the method de- 
tailed in [9] . We compute the corresponding theoretical 
lensing predictions as described in the following section. 
Although these lensing computations follow a fairly stan- 
dard procedure, our procedure is described explicitly be- 
low to make clear which approximations we make. 



B. Lensing calculations 

The lensing measurements that we will use are quoted 
in terms of the variance of the aperture mass M ap , given 
by [57,58]: 
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where J4 is the fourth-order Bessel function of the first 
kind and P K (£) is the convergence power spectrum, given 
by [57-60] 



P K (t) 
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Here P$ l (k, z) is the nonlinear mass power spectrum, and 
a(w) = [1 + z(w)]^ 1 is the cosmic scale factor, w is the 
radial comoving distance defined by 
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and wh = w(oo) is the comoving distance to the horizon. 
The time-variation of the Hubble parameter H is given 
by 



H 



= + zfn m + (1 + z) 2 (i - a„ - O0 + n A . (5) 



fx(w) is the comoving angular diameter distance out to 
w, given by 

K-V 2 smiK-^w) K>0 
f K (w) = { w K = (6) 

(-X)- 1 / 2 sinh((-X)- 1 /2 u; ) K<0 

where K is the spatial curvature defined by 

K = -(^j n k , n k = (1 - n m - n A ). (7) 

The weighting function W(w) is the source-averaged dis- 
tance ratio 



W(w) 



WH y-H, /X fl<(w' — W) , , , . 

Jk(w') 



where G(w) = G(w(z)) is the source redshift distance 
distribution, usually approximated by a fitting function 
of the form 



G{w) 



z T{^) \z 



(9) 



with fitting parameters zq, a and /3. 

The linear power spectrum at redshift z is related to 
the present one by 
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where the linear growth factor g(z) is approximately 
given by [61] 
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Here the parameters O^ and O z K are those at redshift 
z, given by the present values O m and 0\ through the 
relations 

2 

O m (l + z) 3 , (12) 
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H(z)_ 

with H(z) given by equation (5). 

The power spectrum P£ l (k) needed in equation (3) is 
the nonlinear one rather than the linear one Pg(k) given 
by equation (10). Based on a pioneering idea of Hamilton 
et al. [62], a series of approximations [59,63-65] have been 
developed for approximating the former using the latter. 
In terms of the dimensionless power 

4rr 
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the linear power A] on scale k\ is approximately related 
to the nonlinear power A n i on a smaller nonlinear scale 
k u \. We use the Peacock & Dodds' approximation [65], 
where this mapping is given by 

A^(fc n ,)=/»i [A? (A,)] 

and 

k = [l 

with a fitting function 



A 2 .M-X/3 



111 • 



(15) 



(16) 



fnl(x) 



1 + B(3x + (Ax) 



a/3 



1 



(Az)° 3 (0)3 



1//3 



(17) 



parametrized by 

A = 0.482(1 
B = 0.226(1 
a = 3.310(1 
/3 = 0.862(1 
y = 11.55(1 



/o\-0.947 

/o\-1.778 

/q\— 0.224 
Weff/3) 

/o\-0.287 
/o\-0.423 

Here g(0) is the linear growth factor of equation (11) 
evaluated at z = and n c g = d\nPg(k)/d\nk\ is the 
effective logarithmic slope of the linear power spectrum 
evaluated at fci. Since this slope should be evaluated 
for a model without baryonic wiggles, we compute n e s 
using an the Eisenstein & Hu fitting function with baryon 
oscillations turned off. 



(18) 
(19) 
(20) 
(21) 
(22) 



III. EXPERIMENTAL DATA USED 

A. CMB data 

Figure 1 shows the 135 CMB measurements which are 
used in our analysis. Compared to the data set we used 
in [30], we add the new measurements from the Cosmic 
Background Imager (CBI) mosaic [66], the Very Small 
Array (VSA) [67] and Archeops [68]. For CBI, we use the 
year 2000 observations of three pairs of mosaic fields [66] 
but not the deep fields, because it is still unclear whether 
their signal is dominated by CMB or other effects such 
as SZ effect [69]. The Boomerang results updated last 
week [70] and the Acbar results [71] became available too 
recently for inclusion in this analysis, but we do include 
them in the online combined power spectrum described 
below. 
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FIG. 1. CMB data used in our analysis. Error bars do not 
include calibration or beam errors which allow substantial vertical 
shifting and tilting for some experiments (these effects are included 
in our analysis). 

We combine these measurements into a single set of 28 
band powers shown in Figure 2 and Table 1 using the 
method of [30] as improved in [31], including calibration 
and beam uncertainties, which effectively calibrates the 
experiments against each other. Since our compressed 
band powers de are simply linear combinations of the 
original measurements, they can be analyzed ignoring the 
details of how they were constructed, being completely 
characterized by a window matrix W: 



(23) 



where 8Tf = l{l + \)Ci/2tt is the angular power spec- 



trum. 



This matrix is available 



at 
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www .hep.upenn.edu/ ~max / cmb/ cmblsslens .html 
together with the 28 band powers di and their 28 x 28 co- 
variance matrix. The data ^-values and effective ^-ranges 
in Figure 2 and Table 1 correspond to the median, 20th 
and 80th percentile of the window functions W. Com- 
paring Table 1 with the older results from [31], we find 
that the only major change is a shallower rise towards 
the 1st peak due to Archeops, which is able to help cal- 
ibrate Boomerang and other small-scale experiments by 
connecting them with the COBE. Specifically, STg has 
increased by about 10% at I ~ 50 and decreased about 
5% for I ~ 100 — 200, thereby nudging the first peak a 
tad to the right. 
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FIG. 2. Combination of data from Figure 1. These error bars 
include the effects of beam and calibration uncertainties, which 
cause long-range correlations of order 10% over the peaks. In addi- 
tion, points tend to be anti-correlated with their nearest neighbors, 
typically at the level of 10-20%. The curve shows our model best 
fitting CMB+LSS data (second last column in Table 2). 
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B. LSS data 

Measurements of P(k) from Galaxy redshift surveys 
have recently improved in both quality and quantity, and 
the Sloan Digital Sky Survey is set to continue this trend. 
In this paper, we use the power spectrum from the 2dF- 
GRS [72] as measured by [73] . We model the galaxy bias 
as a scale-independent constant 6, and therefore discard 
all 2dF measurements with k > 0.3/i/Mpc to minimize 
our sensitivity to nonlinear clustering and nonlinear bias 
effects. 



Table 1 - Band powers combining the information from CMB data 
from Figure 1. The 1st column gives the £-bins used when com- 
bining the data, and can be ignored when interpreting the results. 
The 2nd column gives the medians and characteristic widths of the 
window functions as detailed in the text. The error bars in the 
3rd column include the effects of calibration and beam uncertainty. 
The full 28 X 28 correlation matrix and 28 X 2000 window matrix are 
available at www.hep.upenn.edu/ ~ max /cmb /cmblsslens. html. 



4 



0.01 0.1 1 

k [l/'h 1 Mpc] 



10 



0[arcmin] 



FIG. 3. 2dF galaxy power spectrum data used. We include 
only points leftward of the dotted line (with k < 0.3h/Mpc) in our 
analysis to stay in the linear regime. The curve shows our model 
best fitting CMB+LSS data (second last column in Table 2). 

C. Lensing data 

For this paper, we use the results from Hoekstra et al. 
[49] , the analysis of 53 square degrees of i?c-band imag- 
ing data from the Red Sequence Cluster Survey (RCS). 
The RCS [45,46] is a 100 square degree galaxy clus- 
ter optical survey designed to provide a large sample 
of optically selected clusters of galaxies with redshifts 
0.1 < z < 1.4 We use only the seven data points which 
have been used in cosmological analysis in [49]. Wc 
model the source redshift distribution, we use the multi- 
redshift part in equation (9), with a — 4.7, /3 = 1.7 and 
zo = 0.302, which are the best fit value given by Hoekstra 
et al. [49]. 



FIG. 4. Lensing data used in our analysis. The two curves show 
our best fit model excluding (bottom) and including (top) lensing 
data. 



IV. PARAMETER CONSTRAINTS 

Our cosmological parameter constraints are summa- 
rized in Table 2. We study three cases: using CMB 
data only, adding LSS information and adding lensing 
measurements too, respectively. These three cases corre- 
spond to the three columns of the table and are discussed 
in the corresponding three subsections below. 
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FIG. 5. Constraints on individual parameters using only CMB 
information. The quoted 1 — a and 2 — a confidence limits are where 
each curve drops below the horizontal dashed lines e _1 /2 ?s 0.6f 
and e~ 2 I' 1 m 0.14, respectively. 

Figure 5 shows our constraints on nine individual cos- 
mological parameters, using only the CMB data. Com- 
paring to numerous previous analyses using CMB data 
alone [8,9,23,30], we see that the addition of new CMB 
data from CBI, Archeops etc. has substantially tightened 
the constraints, especially on the tensor-to-scalar ratio r 
and dark energy J7a- CMB alone is seen to firmly rules 
out models with no dark matter and models with no dark 
energy. (The dark energy constraint depends critically on 
our prior assumption f2t = 0, whereas the dark matter 
result does not.) 

The last four constraints in Table 2 are computed with 
the moments method as in [30] . They are for parameters 
(h, P, Zi on and <jg) that are not independent and merely 
functions of the 11 parameters that define our model grid. 
The Hubble parameter h is given by 



h = 



(24) 



The redshift-space distortion parameter /3 is given by 
[74-76] 



= 



/(On,n A ) 



(25) 



where /(fi m , ^a) ~ ^m 6 i s the velocity-suppression fac- 
tor [75,76] (we evaluate / exactly using the grow A pack- 



age [74] ) . If the diffuse intergalactic hydrogen was reion- 
ized abruptly at a redshift Zi on , then [77]. 



' u\ 2 / 3 

J.9(^ ny>. 



(26) 



in the approximation that z- lon 1 (which we know to 
be the case). The alternative normalization parameter 
erg is the rms fluctuation in an 8/i -1 sphere, given by a 
weighted integral of the matter power spectrum [75] . 



Table 2 - Best fit values and \—a confidence limits on cosmological 
parameters. For the numbers above the horizontal line, the central 
values are the ones maximizing the likelihood (the best fit model). 
For the numbers below the horizontal line, the central values are 
means rather than those for the best fit model [30]. For instance, 
the Hubble parameters for the best fit models are h = 0.68, 0.69 
and 0.61 for the three columns, respectively, which differs from the 
mean values in the table. 
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B. From CMB+LSS data 



C. From CMB, LSS and lensing data 

1. One- Dimensional Constraints 
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FIG. 6. Same as previous figure but combining CMB and LSS 
information. 

Figure 6 shows how the above-mentioned constraints 
tighten up when including the 2dF galaxy clustering 
data. The improvement comes mainly from breaking de- 
generacies involving curvature, dark energy, dark mat- 
ter, baryons and tilt. The f^A-measurement is seen to 
be nicely consistent with the many other pieces of ev- 
idence (SN la etc.) pointing towards O m ks 0.7). The 
2dF bias value b w 0.87 ± 0.06 (interpreted as the real- 
space bias at the effective survey redshift) is consistent 
with that found by the 2dF team [78,79], although on 
the low side (since we are marginalizing over b, this has 
no impact on our constraints on other parameters such 
as as)- The most striking improvement is seen to be for 
the reionization optical depth r, hinting at a detection of 
reionization around z«8. However, as we will discuss at 
length in Section V, the true significance of this detection 
is probably lower than Figure 6 suggests. 



FIG. 7. Same as previous figure but combining CMB, LSS and 
lensing information. 

Above we saw that adding 2dF data to the CMB 
not only tightened parameter constraints, but also pro- 
duced a consistent picture that agreed well with other 
measurements supporting the cosmological concordance 
model. Figure 7 shows that things are quite different 
when adding the RCS lensing data. Although the nomi- 
nal constraints do get substantially sharper as predicted 
[54], most strikingly for the tensor-to-scalar ratio r, the 
harmonious consistency becomes clouded. Adding the 7 
lensing data points raises % 2 by 22.5 (from 38.7 to 61.2 
for about 54 effective degrees of freedom*). Although 
this per se might be considered acceptable, it is striking 
how much the lensing shifts the best-fit values of many 
parameters. Most notably, the reionization optical depth 



"The effective number of degrees of freedom is the number 
of of data points minus the effective number of parameters fit 
for. Although we nominally used 9 parameters, b has an effect 
only when LSS is included and nt has no effect in the absence 
of tensors (and r is indeed consistent with zero) . The effective 
number of degrees of freedom is therefore about 28 — 7 = 21 
for the CMB-only case, 28 + 20 - 8 = 40 for the CMB+LSS 
case and 28 + 20 + 7-8 = 47 for the CMB+LSS+lensing case. 
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t increases from 0.06±0.03 to O^O+O.OS 1 and the best fit 
cosmological constant Ha drops from 0.71 to 0.55, which 
sits uncomfortably with, e.g., SN la constraints. The cos- 
mological model is pushed into an uncomfortably config- 
uration by the lensing data, creaking and squeaking but 
refusing to outright collapse. 

As we will see in Section V, this tension comes from 
the overall normalization of the power spectrum, with the 
RCS lensing data pushing the best-fit cr 8 value up from 
0.71 to 0.91. This normalization agrees with that of some 
other cosmic shear surveys [47,49-51], but not with that 
preferred by CMB and LSS data. In Section V, we will 
comment on uncertainties and other lensing data with 
different normalization. 



2. Two- Dimensional Constraints 

For more detailed insight into our parameter con- 
straints and degeneracies, we now turn to joint con- 
straints on pairs of parameters. 



CMB alone 



FIG. 8. Constraints in the (r, r)-plane. The large pink/light 
grey region is ruled out by CMB alone at the 95% confidence level. 
The red/grey area shows the region excluded by adding additional 
constraints from 2dF. The black curve shows the allowed region 
when combine CMB, 2dF and RCS lensing. 



^The Gunn-Peterson constraint implies that our Universe is 
highly ionized for z < 6, and the recent detection of a Gunn- 
Peterson trough in z ~ 6 quasars [80] has been interpreted 
as the tail end of reionization. Moreover, it has been argued 
that the temperature of the intergalactic medium at z ~ 3.4 
inferred from He II reionization requires z < 9 [81,82], nicely 
consistent with our best fit CMB+LSS value z ps 8. How- 
ever, models will less abrupt reionization still allow substan- 
tially larger r-values consistent with our CMB+LSS+RCS 
constraint r = 0.20 ± 0.03 [33,83,84], so the issue of whether 
to trust the lensing normalization cannot yet be settled by 
reionization constraints. 




0.4 0.6 



FIG. 9. Constraints in the (n s ,r) plane. The large pink/light 
grey region is ruled out by CMB alone. The red/grey region be- 
comes ruled out when adding 2dF information. The black curve 
shows the allowed region when combine CMB, 2dF and RCS lens- 
ing. 

Figure 8 shows that 2dF breaks the CMB degeneracy 
between reionization and gravity waves and produces es- 
sentially independent (separable) constraints on r and r. 
The gravity wave limits agree well with those found by 
[36] who also included 2dF and Archeops data. The fact 
that the allowed regions excluding and including lensing 
data have essentially no overlap illustrates the poor con- 
sistency between CMB+LSS and lensing data sets men- 
tioned above. 

Figure 9 shows constraints in the parameter space 
(n s ,r) where inflationary models can be tested. These 
limits from CMB only and from CMB+LSS are compa- 
rable with those in older work requiring stronger priors. 
For instance, comparison with a corresponding plot from 
a year ago [30] shows new data has squeezed the con- 
straints substantially, particularly on the gravity wave 
amplitude. * Now CMB data alone gives r < 0.48 at 95% 
level, obtainable only be combining CMB and LSS data 



*As discussed in [85], there has been a fair amount of nota- 
tional confusion in the literature surrounding the tensor-to- 
scalar ratio r. There are two logical ways to define this ratio: 
either in terms of the fundamental parameters of the power 
spectrum (or, equivalently, of the inflationary model space), 
or in terms of the observables, usually the CMB quadrupoles. 
As in [30], we adopt the former approach and define 



A s 



(27) 



where A s and At are the scalar and tensor fluctuation am- 
plitudes as denned in [85]. For inflation models where the 
slow-roll approximation is valid, this ratio is related to the 
tensor tilt nt by the so-called inflationary consistency condi- 
tion [85,86] 
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back then. Now CMB+2dF lowers the limit to r < 0.34. 
The lensing data nominally helps even helps more, tight- 
ening the constraint to r < 0.13 — we return to the issue 
of whether this is believable in Section V. 
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FIG. 10. Constraints in the (u,j m! W] 3 )-plane. The large 
pink/grey region is ruled out by CMB data alone at 95% confidence. 
The smaller red/shaded region becomes excluded when imposing 
additional constraints from 2dF. The thick black curve shows the 
allowed region when combine CMB, 2dF and RCS. The rectangu- 
lar band by the thin straight lines shows the BBN constraints from 
[87]. 



with CMB and galaxy clustering data, we have perform 
a detailed 9-parameter analysis of recent lensing, CMB 
and galaxy clustering data. 

For the CMB part, we constructed and used a data set 
combining all available data taking calibration errors and 
beam uncertainties into account, effectively calibrating 
them off of each other. Using this combined power spec- 
trum for a "MAP versus world" comparison next month 
will provide a powerful test of how realistic the error es- 
timates have been in the CMB community. 

Even though 2002 has added substantial amounts of in- 
formation on the CMB power spectrum, e.g., from CBI 
and Archeops, the combined CMB and 2dF data remains 
comfortably consistent with a simple flat adiabatic scale- 
invariant concordance model with Cl\ = 0.72 ± 0.09, 
/i 2 O cdm = 0.115 ± 0.013 and h 2 n h = 0.024 ± 0.003. The 
new data has added a hint of reionization around z ~ 8. 
These constraints assumed negligible contributions from 
spatial curvature and massive neutrinos. 

Including the RCS lensing data in our analysis con- 
firmed with real data what was predicted by Fisher ma- 
trix forecasts [54] : that lensing has the statistical power 
to substantially help with cosmological parameter con- 
straints in two ways: by breaking degeneracies and by 
providing cross-checks. Intriguingly, one of the cross- 
checks failed, so let us now turn to this failure in more 
detail. 



Let us now turn to the cosmic matter budget. A year 
ago [30,23], the CMB+LSS constraints on the baryon 
density Wb were beautifully consistent with those from 
Big Bang Nucleosynthesis (BBN), u h = 0.02 ±0.002. Al- 
though the agreement is still acceptable, Figure 10 shows 
that the CMB+LSS baryon density is creeping back up 
again towards the higher range that was favored two 
years ago [1,2]. This slight tension is reflected in x 2 in- 
creasing by 3 if we add a the BBN prior to our likelihood 
analysis. 



V. DISCUSSION 

Since cosmic shear measurements have now improved 
to the point where they deserve to be treated on par 



200 



-n t . 



(28) 



A common alternative definition of the tensor-to-scalar ratio 
is the quadrupole ratio 



R ■: 



^scalar ' 



(29) 



Writing the ratio r/R is typically between 2 and 5 — it de- 
pends on the values of Qa and fi^ via the late integrated 
Sachs- Wolfe effect. 



A. The (78 problem 

Our results suggest that the tension between the lens- 
ing and other data comes from the overall normalization 
of the power spectrum. We saw that the CMB+2dF data 
preferred as = 0.74, in good agreement with results from 
last year [78]. In contrast, the RCS lensing data prefers 
higher values [49] and pushed the normalization up to 
as = 0.92 in our joint analysis. 

Physically, it is easy to see why this pushed other pa- 
rameters in the way it did, notably increasing r and de- 
creasing 57a- If we increase the value of as, the CMB 
and galaxy power spectra will both increase in amplitude. 
This causes no problems for fitting the 2dF data, since 
our method can merely readjust their bias, but causes se- 
rious difficulties with the CMB. The acoustic peaks can 
be brought back down to acceptable heights by increasing 
t, but this provides only a partial solution to the prob- 
lem, since it leaves the power on very large I < 10 scales 
larger than observed by COBE DMR. Lowering £1a helps 
with this, since it increases the linear growth factor and 
thereby lowers the entire CMB power spectrum for fixed 
current matter clustering. Q\ can only be lowered by 
a moderate amount, however, since the matter density 
~ 1 — 57a will increase and give the galaxy power 
spectrum the wrong shape. The response of the parame- 
ter fitting algorithm is therefore to do a little bit of both 
tweaks, increasing r and lowering Indeed, similar 
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results were obtained in a smalled parameter space by 
Hoekstra et al. [49] , who found that they could only rec- 
oncile their lensing normalization with CMB for a rather 
large reionization optical depth r = 0.12 ± 0.04. 

To verify that this is the correct physical interpreta- 
tion of what is going on, we repeated our entire analysis 
with the amplitude of the RCS data lowered by a fac- 
tor of (0.7/0.86) 2 , which is roughly the ratio of normal- 
izations preferred by CMB+2dF and RCS respectively. 
This analysis gave results fully consistent with those 
from CMB+2dF alone, for example, r = 0.05 ± 0.03, 
Ov = 0.76 ±0.13, u dm = 0.11 ±0.02, uj h = 0.024 ±0.005, 
n s = 1.00 ± 0.04, A s = 0.52 ± 0.10, r = 0.10 + 0.23, 
b = 0.82 ± 0.09, z lon = 7.64 ± 3.18 and a 8 = 0.70 ± 0.11. 
These results not only demonstrate consistency, but also 
show that by combining lensing and CMB and LSS, we 
can get stronger constraints on parameters like ft a, n Sl 
as and Zi on , just as had been predicted using information 
theory [54]. 

Table 3 - Recent measurements of the power spectrum normal- 
ization as- For results quoted as fitting formulas involving Q 1 1 1 , we 
have set Q m = 0.3. Error bars are la — for [49] and [52], we have 
simply divided their quoted 95% errors by 2. 



Analysis 




°-8 


Clusters: 






Pierpaoli et al. (2001) 


[88] 


1 02+ ' 07 


Borgani et al. (2001) 


[89] 


„ _„ + (). 1)8 
°' 76 -0.05 


Reiprich & Bohringcr (2001) 


[90] 


n fis+ ' 08 

U - bS -0.06 


Seljak et al. (2001) 
Viana et al. (2001) 


[91] 
[92] 


0.75 ±0.06 
0.61 ± 0.05 


Bahcall et al. (2002) 


[93] 


0.72 ±0.06 


Pierpaoli et al. (2002) 


[94] 


n 77 +0.05 


Weak lensing: 




n 71+O.O6 


Jarvis et al. (2002) 


[52] 


Brown et al. (2002) 


[95] 


0.74 ± 0.09 


Hoekstra et al. (2002) 


[49] 


n »R+ ' 04 


Van Waerbekeei al. (2002) 


[51] 


0.97 ± 0.06 


Baconet al. (2002) 
Refregieret al. (2002) 


[47] 
[50] 


0.97 ±0.13 
0.94 ±0.14 


CMB: 






Lahav et al. 2001 


[78] 


0.73 ± 0.05 


Mclchiorri & Silk 2002 


[3d] 


0.70 ±0.05 


Lewis et al. 2002 


[32] 


0.79 ± 0.06 


This work 




0.74 ± 0.06 



So is the CMB normalization too low or is the lens- 
ing normalization too high? As shown in Table 3, as is 
emerging as the currently most controversial cosmological 
parameter, and it will be crucial to use further cosmolog- 
ical data to get to the bottom of this issue. The table 
illustrates what has been emphasized by many authors, 
namely that the scatter between as measurements using 
different techniques is substantially larger than the for- 
mal error bars, triggering unpleasant flashbacks to the bi- 
modal and bitter debate about the Hubble parameter h. 
Although much work clearly needs to be done to resolve 
this issue, there are some encouraging indications that 
the as gap may be closing. The CMB normalization has 
increased slightly in the last year, with the first acoustic 



peak growing by 10% between [30] and our present com- 
pilation. CMB experiments can now be more effectively 
calibrated off of each other, and Archeops connects the 
calibration of DMR to that of the many experiments sen- 
sitive to the 1st acoustic peak and beyond. The cluster 
normalization has dropped somewhat because of recali- 
brations of the mass-temperature relation (see Table 3). 

Although the lensing normalization has been uniformly 
high across groups [47,49-51], the largest and most recent 
data set released has now produced the much lower nor- 
malization as ~ 0.71 [52] (for f2 m = 0.3). Indeed, the 
Hoekstra et al. normalization may also be coming down 
by 5-10% (Hoekstra, private communication), partly be- 
cause of a new and more accurate version of the Peacock 
& Dodds fitting formula [65]. 

Moreover, our analysis did not marginalize over the 
source redshift parameters a = 4.7, (3 = 1.7 and zq — 
0.302 from equation (9), but merely used the best fit 
values from [49] . Such marginalization would be expected 
to somewhat weaken the as constraints that we obtained 
from the RCS lensing data. As an extreme example, we 
recomputed the shear power spectrum P K (£) replacing 
the characteristic source redshift parameter z s = 0.302 
with the 3 — a lower and upped limits on its value from 
[49]. These extreme values z s = z = 0.274 and z s = 0.337 
made the shear fluctuation amplitude 10 — 15% lower 
and higher, respectively, so lowering the source redshifts 
at the 2a level would close approximately half of the as 
gap between the RCS and the CMB. Despite the current 
discord, there is thus real hope that a beautiful consistent 
picture of cosmology will eventually emerge. 

B. Are the constraints too good to be true? 

In addition to the above-mentioned tension regarding 
the best fit values of some parameters, there is a slight 
puzzle regarding their error bars: some constraints seem 
a bit too good to be true. Specifically, the CMB er- 
ror bars in Table 2 arc in some cases comparable to 
those forecast for MAP [96] using the Fisher matrix 
technique [97]. This puzzling fact is not limited to the 
present paper, but appears rather generic. For instance, 
consider the scalar spectral index n s . The forecast for 
MAP (2 years of data including polarization) is accuracy 
An s = 0.11 with no prior assumptions, An s = 0.045 
assuming Ok = and An s = 0.021 when adding SDSS 
galaxy clusterinng data (the final luminous red galaxy 
sample) [96]. We found An s ~ 0.06 from current CMB 
alone (Table 2, assuming Qk = 0), comparable to other 
studies in the recent literature. For instance, Lewis & 
Bridle [32] obtain An s = 0.04 including 2dF data, and 
the Boomerang team reports An s = 0.06 — 0.08 for vari- 
ous priors [70]. 

To gain insight into this puzzle, we compute the Fisher 
matrix [97] corresponding to the current data. In the 
approximation that all the CMB information is coming 
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from the mean rather than the covariance of the power, 
the equations in [97] give 

F cmb = D'WS^W'D, (30) 

where £ is the 28 x 28 covariance matrix of our combined 
CMB power spectrum measurements from Section III A, 
W is the corresponding £ max x 28 window matrix and the 
f mal x 28 matrix 

D« = ±Slf (31) 

gives the derivatives of the power spectrum with respect 
to the cosmological parameters pi. We compute an anal- 
ogous Fisher matrix F 2c if using the covariance matrix 
and window functions for the 2dF power spectrum from 
[73]. We then evaluate the attainable error bars from 
CMB alone as A Pl = {F'^fJ 2 and from CMB+2dF 

as Api = ([F cm b + F2df] -1 )i/ 2 - For * ne sca l ar spectral 
index with an = prior, this gives An s = 0.16 from 
CMB alone and An s = 0.13 for from CMB+2dF, a factor 
2-3 larger than the above-mentioned likelihood analysis 
results from us and others. 

Comparing with our Fisher matrix results, the most 
suspicious error bar in Table 2 is that on the optical depth 
r, where we quoted At = 0.06 from CMB alone and the 
Fisher matrix gives At = 0.18. For the t case, there 
is wide scatter in the literature. On the optimistic side, 
the CBI team reported r = 0.09t 1 Q7 [29], and a strong 
detection of reionization was claimed from earlier data 
[98]. On the pessimistic side, the 2dF and Boomerang 
teams report a 2a limit r < 0.5 [23,70], in line with the 
Fisher forecast. 

If our t errors are indeed too small due to some artifact 
of our method, this would weaken the constraints on the 
normalization as, since CMB constrains mainly the com- 
bination a$e~ T . However, the conclusion that er 8 = 0.9 
forces the uncomfortably high value t w 0.2 would still 
stand, since it follows directly from this scaling relation. 

Efstathiou & Bond [99] have shown that there are 
physically understandable reasons why the Fisher ma- 
trix approach can significantly overestimate the errors on 
cosmological parameters, especially when strong degen- 
eracies are present and that the curvature of a banana- 
shaped allowed region becomes important. This is indeed 
the case with t, with its strong degeneracy with the nor- 
malization parameter A s , and a rigorous Bayesean confi- 
dence interval taking into account the asymmetry t > 
(neglected in the Fisher treatment) also matters. "Weak" 
priors can also have a strong effect, notably on h and 
gravity waves. Hopefully, our concerns will therefore go 
away when the MAP data arrives, in particular when all 
major degeneracies are broken with other datasets and 
priors. Since equation (30) is trivial to apply, it would 
be prudent to include a Fisher analysis of the data used in 
all future parameter constraint papers, as a reality check. 



C. Outlook 

In conclusion, the pace at which new measurements 
have tightened up constraints on cosmological parame- 
ters is breathtaking. The constraints become particularly 
tight if one uses theoretical prejudice to impose spatial 
flatless Ok = like we have done, interpreting the CMB 
indication « as support for inflation and perfect 
flatness and thereby breaking the worst parameter de- 
generacy of all by fiat. It can be argued that precision 
cosmology is already here, before MAP, and now faces 
its first baptism of fire: will the current precision mea- 
surements of cosmological parameters agree with MAP? 
Next month we will know whether the whole edifice goes 
down in flames or stands stronger than ever. 
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